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Magnetization plateaus in the Shastry-Sutherland model for SrCu2(BOs)2 are studied by the 
perturbation expansion method. The fourth-order effective Hamiltonian which describes the 
dynamics of triplet dimers (TD's) with S* ot = 1 in the singlet sea is derived and then partially 
diagonalized for the space that consists of the TD configurations with the lowest second-order 
energy. The fourth-order terms are treated within a low-density approximation. Our procedure 
makes clear how TD interactions are responsible for the formation of magnetization plateaus. 
Particularly, the 1/4-plateau is obtained by the fourth-neighbor TD repulsion in the fourth-order 
perturbation, and a diagonal stripe arrangement of TD's appears at this plateau. 
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§1. Introduction 

Two-dimensional frustrated Heisenberg antiferromag- 
nets have been considerable interests for a long time be- 
cause of unique magnetic properties which reflect un- 
usual low-energy spectra due to a strong geometrical 
frustration. A typical example of such frustrated an- 
tiferromagnets is the spin-1/2 Kagome antiferromagnet, 
where a disordered spin-liquid like ground state, a contin- 
uum of singlet excitations adjacent to the ground state, 
and gapped triplet excitations are obtained by numerical 
investigations & 

For the last two years, there has been a growing in- 
terest to a two-dimensional frustrated antiferromagnjet 
which is called the Shastry-Sutherland model (SSM)J$ 
because Kageyame et al. and Miyahara and Ueda found 
that the SSM is realized in SrCu2(B0 3 ) 2 BS The Hamil- 
tonian of the SSM under the uniform magnetic field H 
is defined by 

(iJ) «i>J» 

where the first (second) term represents the intradimcr 
(interdimer) coupling, and M — J2i$i- (See Fig- 1.) 
Strength of the interdirnex coupling in SrCu 2 (B03)2 is 
obtained to be A ~ 0.63.0'B' A distinctive feature of this 
model is that the orthogonal dimer structure shown in 
Fig. 1 leads to the exact dimer singlet ground state for 
A < A c , and_prohibits one-TD hopping up to the fifth 
order in AJjlP As for the quantum phase transition at 
A = A c , Koga and Kawakami obtained A c = 0.68 and 
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found that the exact dimer singlet ground state becomes 
to be unstable against the plaquette RVB state there.EP 




Fig. 1. Structure of a two-dimensional network formed by Cu + 
(closed circles) in SrCu2(B03)2- The vector p\ (<5;) is the unit 
vector parallel (perpendicular) to the axis of a dimer at site I. 
The arrows indicate first-neighbor (In), second- neighbor (2n), 
third-neighbor (3n, 3n') and fourth-neighbor (4n) dimer pairs. 
Note that there exist two types of third-neighbor dimer pairs. 
For brevity, we call that indicated by "3n" the third- neighbor 
pair. 



Low-energy excitation spectra in the SSM have been 
investigated extensively because it is expected that the 
orthogonal dimer structure makes the spectra quite 
unique: the localized nature of one TD leads to a nar- 
row band of scattering states of two TD's, while a TD 
pair with a short interdimer distance should have inher- 
ent dynamics in the multiplicityocP For instance, first- 
neighbor TD pair excitations can propagate through the 
third-order perturbation. As for experimental studies, in 
addition to the spin gap of 3.0 meV, the second spin gap 
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of 4.7 meV was observed hy ESR and neutron inelastic 
scattering experiments S3 Elf In particular, the neutron 
inelastic scattering experiment revealed that the second 
lowest triplet branch shows dispersive behavior although 
the lowest triplet branch is almost flat. In a Raman 
scattering experiment, a singlet excitation mode at 3.7 
meV was observed^ These singlet and triplet bound 
states are considered to originate from first-neighbor TD 
pair excitations. On the other hand, first-neighbor TD 
pair excitations in the quintet sector lead to antibound 
states, and the lowest quintet branch consists of second- 
neighbor TD pair excitations. The lowest quintet branch 
shows rather weak dispersion, and the binding energy is 
much smaller than those of the singlet and triplet bound 
states. 

The orthogonal dimer structure also gives rise 
to an unique feature of the magnetization pro- 
cess: SrCu2 (1303)2 shows intermediate magnetization 
plateaus at M = Af s /4 and where M s repre- 

sents the saturated magnetization^ To understand these 
plateaus, several groups derived the third-order effec- 
tive Hamiltonian describing the dynamics of TD with 
gtot = 1 [ n the singletsea. and analyzed the resultant ef- 
fective Hamiltonian.EJ'EJ'EjE^P These investigations pre- 
dicted the appearance of a l/JLplateau, which was ob- 
served in a recent experiment £3 However, the origin of 
the 1/4- and 1/8-plateaus have not been clarified yet. 
Three scenarios to understand these two plateaus have 
been suggested. The first one is that these plateaus are 
produced by long-range TD repulsions due to higher- 
order perturbation. The second one is that some spin 
interactions omitted in the SSM play an essential role to 
realize these plateaus. For instance, it was shown that 
a magnetization plateau at M — M s /A is obtained in a 
generalized SSM.EZP The third one is that the relevant 
cluster to understand these plateaus is not a spin dimer, 
but a spin plaquette, which is based on the fact that 
SrCu2(BQs)2 locates near the dimer-to-plaquette phase 
boundary.© 

In the present paper, we show that the fourth-order 
perturbation explains the 1 /4-plateau invoking no addi- 
tional spin interactions. We also describe in detail our 
previous brief report in Ref. 15 concerned with the third- 
order perturbation. 

The present paper is organized as follows: in § 2, we 
derive the fourth-order effective Hamiltonian with a low- 
density approximation on the basis of the coupled dimer 
picture of the SSM. The fourth-order effective Hamil- 
tonian is partially diagonalized for the subspace with 
the lowest second-order energy, and then the resultant 
Hamiltonian is analyzed by the numerical diagonaliza- 
tion method. In § 3, we give discussions about our cal- 
culation. It is particularly discussed how lower-order flip 
terms affect higher-order TD repulsions. In § 4, we sum- 
marize our results. 

§2. Perturbation Expansion 



2.1 Zeroth- Order Perturbation 

We start with the dimer limit. Eigenstates of a dimer 
are the singlet state \s) and the triplet states \t\), |io)> 
|£-i), where the suffix denotes the total S z . The lowest 
energy states with magnetization M(> 0) for the whole 
system consist of |s) and \t\). We assume that the low- 
est energy states are in the degenerate space excluding 
\t ) and when the interdimer coupling is taken into 
account x3 Then we calculate the fourth-order effective 
Hamiltonian 



-1/4PE 



where 



M(A 



Err 



4PE 

sg 



H) 



3iV D 



(2) 



(3) 



is the absolute ground state energy of the system with 
Nu dimers, and 



A 3 



A J 



At PE = 1 - A 2 - ^ - ^ (4) 
sg 2 8 w 

is the fourth-order series of the spin gap. In eq. (||), 
V n represents interaction among TD's in the nth-order 
perturbation. If the TD interaction leads to a kink in 
the minimum interaction energy per dimer, e(m), as a 
function of density of TD's, m{= M/M s ), then a mag- 
netization plateau is realized. 

To write down the explicit form of V n , it is convenient 
to introduce a pseudo-spin operator Ii for each dimer 
site. The up and down spin states of Ii correspond to 
the original dimer states as follows: 



IT>! "!*!>!. li>J 



(5) 



We also define the number operator of up spin and the 
spin-flip operator: 



■IT if 



(0) 



(7) 



In terms of the original dimer states, ui means the num- 
ber operator of TD and fi t i> the hopping operator of TD. 

2.2 First-Order Perturbation 

The first-order interaction, Vi, is calculated as follows: 



(M'>£ln 



umv 



(8) 



which is the nearest-neighbor Ising model on a square 
lattice. As is well known, TD's can be arranged on the 
square lattice avoiding first-neighbor TD pairs for m < 
1/2, but can not for m > 1/2. The first-neighbor TD 
repulsion therefore leads to a kink in e(m) at m = 1/2 
as shown in Fig. 2, which means the appearance of 1/2- 
plateau. The Neel type arrangement of TD's is stabilized 
at the 1/2-plateau. The lower critical field of the 1/2- 
plateau is given by A^ E (= 1). The most important fact 



-F2n— >ln — FJ 



(13) 



and 



Fln<->ln — ^ Uifi-S u l+8i, 



(14) 



Fig. 2. First-order approximation of the minimum interaction en- 
ergy eipe("2) and TD configuration in the Neel state. The 
dimers with open (closed) circles are in the triplet (singlet) dimer 
state. 



We are going to use the same argument to study if 
higher-order TD repulsions lead to additional plateaus. 
The situation becomes more and more complicated as 
the order of the perturbation becomes higher: in H^^, 
there exist many types of diagonal interactions and off- 
diagonal interactions. To overcome this difficulty, we 
utilize the fact of V± being a diagonal operator, which 
enables us to carry out partial diagonalization of H^ E . 
To be precise, we direct our attention to the subspace 
without first-neighbor TD pairs for m < 1/2, in which 
the first-order energy is the lowest value. We write the 
new effective Hamiltonian as follows: 
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HiP = E e + M(AgJ E - H) + J2 * n V n for M < ^, 

n=2 

(9) 

where the constraint excluding first-neighbor TD pairs 



uiUi± Pl = for any / 



(10) 



is imposed. The bar on operators means that the oper- 
ators are defined in the constrained space. 

2.3 Second- Order Perturbation 

The second-order interaction, V2, is calculated as fol- 
lows: 

y 2 = \ um ' + \ 51 uiuv ~ 9 Y u i u i+pi u i-pi 



(M')eln (Z,i')e3n I 

-|(-Fln— >2n + i"2n— >ln + -Fln^ln)- 



(11) 



Here, 



are correlated hopping terms: i*i n — >2n (-P211— »in) annihi- 
lates a first (second) neighbor TD pair and creates a sec- 
ond (first) neighbor TD pair, and Fi n ^i n annihilates a 
first-neighbor TD pair and creates another first-neighbor 
TD pair. As we will see below, i*i n _>2n and -Fan-tin 
play a crucial role to produce an intermediate plateau 
at m = 1/4. 

As mentioned previously, we calculate V2 in the sub- 
space without first- neighbor TD pairs. Noting that the 
matrix elements of the three-body and correlated hop- 
ping terms are zero if initial and/or final states contain 
no first-neighbor TD pairs, we get 



1 



V 2 = - 



2 ^ 

(l,l')e3n 



UlUl'. 



(15) 



This equation shows that the third-neighbor TD repul- 
sion dominates the low-energy physics in the second- 
order perturbation. 

The third-neighbor TD repulsion in V2 increases the 
energy of the Neel state, which is the lowest energy state 
at m = 1/2, by A 2 7Vd/4. Thus the second-order pertur- 
bation leads to a new kink in e(m). To determine the po- 
sition of the kink, we calculate the second-order approx- 
imation of the minimum interaction energy, £2pe(' 7 i)j for 
various finite size clusters up to iV D = 36. The shapes 
of clusters used in the present work are given in Fig. 3, 
where the periodic boundary condition is imposed. We 
show the result of £2Pe(w) in Fig. 4. It is found that the 
kink appears at m = 1/3. The lowest-energy state at 
m — 1/3 is uniquely determined as the 1/3-stripe state 
which is defined in Fig. 4. This result is the same as the 
results obtained by Miyahara and Uedat3 and Momoi 
and Totsuka.BEJ 



Equation (|15|) shows that no off-diagonal terms exist 
in the second-order interaction, which means that higher 
order terms may play an important role. Thus we can 
carry out partial diagonalization again. Then we study 
the subspace without first- and third-neighbor TD pairs 
for m < 1/3, which is the lowest energy states within 
the second-order approximation. We write this effective 
Hamiltonian in the following form: 



o}4PE 771 

n cft - tig 



M(Agf E - H 
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^A"K for M < ^f-, 

n — 3 

(16) 



with the constraint 



UlUl± pi = Ul +pi Ul- f 



for any I. 



(17) 



The tilde on operators indicates the exclusion of first- 
and third-neighbor TD pairs. 

As long as the interdimer coupling A is small enough, 



ory. This gives a support to the validity of our approach 
to SrCu 2 (B0 3 ) 2 . 




Fig. 5. Exact diagonalization results of the minimum interaction 
energy, £ed(™)i for the 12-dimer cluster. 



Fig. 4. Second-order results of minimum interaction energy, 
e 2PE( m )i at A = 0.63, and TD configuration in the 1/3-stripc 
state. 



the lowest energy states for to < 1/3 are in the sub- 
space without first- and third-neighbor TD pairs due to 
TD repulsions up to the second-order perturbation. On 
the other hand, when A is increased, the energy-gain due 
to TD hopping becomes to be more important than the 
TD repulsions. There should be an energy level cross- 
ing between these two regions. Our assumption for the 
restriction of the subspace is expected to be valid for 
wider range of A as the density of TD's is lower, because 
the TD hopping in low-order perturbation is achieved by 
the correlated hopping. To check the validity of our per- 
turbation theory for SrCu2(B03)2 with A ~ 0.63, we use 
the exact diagonalization method and calculate the mini- 
mum interaction energy, eed(w), on the 12-dimer cluster 
which matches the structure of the Neel and 1/3-stripe 
states. The result is shown in Fig. 5. We find no en- 
ergy level crossing up to A ~ 0.63, which means that the 
lowest energy states in SrCu2 (1303)2 are adiabatically 
connected to the states retained in our perturbation the- 



2.4 Third-Order Perturbation 

We turn to the third-order interaction, V3. The cal- 
culation of this term is somewhat tedious, but can be 
achieved straightforwardly. The result is given by 

v 3 =-- u i u i' + 4 zl u i u i' + 4 E um ' 

(l.l')eln (Z,2')62n (Z,Z')e3n 

o 4 o 

+ (three-body terms). (18) 

We find in V3 that second-neighbor TD repulsion appears 
for the first time in the third order. We also note that 
the coefficient of Fi n _ >2 n is different from that of F 2n ~^ in, 
i.e., V3 is not an Hermite operator. All the off-diagonal 
terms in V3 take the form of correlated hopping. The 
explicit expressions of the three-body terms are given in 
Ref. 16. The three-body terms, however, play no role 
in our treatment, because there are no matrix elements 
of the three-body terms in the subspace without first- 
neighbor TD pairs and we will use the two-body approx- 
imation in the fourth-order perturbation. 
Using eqs. (||), (11) and (18), we obtain 

^3 = ^J2(ui+ Pl + Ul +Pl ){u l+Sl Ui-8 l +fl +Su l-6 l )- (19) 

Note that the second- neighbor repulsion exists in V3, but 
does not in %. The processes which cancel the second- 
neighbor TD repulsion are shown in Fig. 6. We assume 
there exist two TD's at l+pi and l+5i in the initial state. 
Operating (A 2 /4)i*2iwin on the initial state, we get the 
intermediate states with a first-neighbor TD pair, whose 
energy is higher by AE — A/2 than the initial state. Op- 
erating (A 2 /4)Fi n ^2n on the intermediate states, we ob- 
tain the final state which is the same as the initial state. 
These second-order processes give a second-neighbor TD 
attraction — (X 3 /A)ui+ Pl ui + s, , which cancels the second- 
neighbor TD repulsion in the first-order process. 
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Fig. 6. Second-order processes which cancel the second-neighbor 
TD repulsion in V3. The number attached on each arrow repre- 
sents the matrix clement. 



The third-order diagonal term in V3 , which is a three- 
body repulsion, does not increase the energy of the 1/3- 
stripe state. Thus we obtain no plateau originating in 
the third-order effect. From this order, there exist cor- 
related hopping terms. To study the effects of the cor- 
related hopping terms on the magnetization curve, we 
diagonalizc numerically V3 on finite size clusters with 
Nn = 24, 36. The resultant magnetization curve at 
A = 0.63 is shown in Fig. 7 together with the exact 
magnetization curve for the cluster with 7V D = 12. The 
correlated hopping terms result in a finite width of the 
magnetic field region between m — and 1/3. However, 
this region is much narrower than that in the exact re- 
sult. This result indicates that the energy-gain due to 
the third-order correlated hopping is rather small. So 
we expect that a magnetization plateau appears if there 
exists a TD repulsion responsible for a magnetization 
plateau in the next order. 

0.35 
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- 3PE (JV D =24) 
3PE (iV D =36) 
■ ED (iV D =12) 



Fig. 7. Magnetization curve at A = 0.63. The solid and long 
dashed lines are the results of the third-order perturbation 
(3PE), and the short dashed line is the result of the exact diag- 
onalization (ED) of the SSM for the 12-dimcr cluster. 



2.5 Fourth- Order Perturbation 

Up to now, all the perturbation processes have been 
taken into account in our calculation. In the fourth-order 
perturbation, such calculation is far from transparent. 
So we invoke the fact that we are now interested in low 



density region of TD's, m < 1/3, and neglect three or 
more body terms in V4. 

The fourth-order term V4 takes the following form: 

(M'}61n (;,Z')S2n (M')G3n 

7^-Fliw2n + ^-F2n_>ln + --Fln<->ln 



UlUp 

'(M')e4n 



1 \ - 

+ 8 £ 

<U'}e 

+ JqY^[( u I+2Pi+Si +Ul+S l -Ul +2 p,-5, -«/-«,) 

xfi,i+ Pl - (pi -» -pi)] 

Xfl,l+pi+5i + (Si -» Si)] 

+ J^^2[ui+p l (.fl,l+2p l +S l - fl,l+2 Pl -Si) - (pl —> -Pi)] 
I 

~~l^2( I l~ I l^-p l +6, I l+p, I l+5, + n ' C ') 
/ 

+ 32 

,l+pi+Sifl+pi,l+Si 

+(three or more body terms). (20) 

We find in V4 that fourth-neighbor TD repulsion appears 
for the first time in the fourth order. The matrix element 
of the second-neighbor TD repulsion is three times as 
strong as that of the fourth-neighbor repulsion. As for 
the off-diagonal part, we find that TD pair hopping is 
possible from this order. 

Calculation of matrix elements in the two-body ap- 
proximation of V4 is easily carried out by introducing 
the two-TD basis set in the momentum space, because it 
is reduced to partial diagonalization of a finite size ma- 
trix due to luck of the one-TD hopping. The result is 
given by 

Vj - g ^ uiu v 

(M')G4n 



-^2,[f 1,1+5^1+ p l -S l 



Ul- 



- Sl )-(Si 



-Si) 



+ 16 

,i+2p l (ui+ Pl +s l +ui +Pl - Sl ) 



~IQ^ fl,l+25i (h+pi+Sj +^l-p,+S l ) 

+ 16 ^2ifl,l+2pi+26iUl +pi+ 26i + (Si —> 
I 

^J^Z^^H-Pi+sJl+pJl+8, + n ' c ')- 



Si)] 



(21) 



Note that the second-neighbor TD repulsion exists in 
V4, but does not in V4. The processes which cancel 



Fig. 9. Fourth-order results of minimum interaction energy 
e^ pE (m), for the diagonal and low-density approximation of 
W*j B , and TD configuration in the 1/4-stripe state, where we 
use A = 0.63. 



1 / 4-plateau observed experimentally. The width of the 
1/4-plateau is also consistent with that of the experi- 
ment. 
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Fig. 10. Magnetization curve at A = 0.63 calculated by using the 
perturbation expansion method up to the second order (2PE), 
the third order (3PE) and the fourth order (4PE). Calculations 
of the third- and fourth-order perturbation are carried out for the 
i/32 X vTS cluster (Ny> =24). Tba circles represent experimental 
data obtained by Onizuka et a/.JU where we use the intradimer 
coupling 84 K and the g-value 2.05. 



§3. Discussions 

In our perturbation theory, we first calculate 7i*g E , 
and then calculate 'H*g E for m < 1/3 in the subspace 
without first- and third-neighbor TD pairs. We here 
comment on the TD repulsions in these two Hamilto- 
nians. In Fig. 11(a), we show the two-body repulsions 
in T~Ccff E - There exist first-, third-, second- and fourth- 
neighbor TD repulsions in 7i^f E , which appear in the 
first-, second-, third- and fourth-order perturbation for 
the first time, respectively. At X — 0.63, the matrix el- 
ements of these repulsions are given by 0.39, 0.41, 0.12 
and 0.02 for first-, third-, second- and fourth-neighbor 



Fig. 12. Interaction energy, e'(m), and the lowest energy config- 
urations at m = 2/9 and 1/6 for "H 1 ^, where we use A = 0.63. 
The correspondence of the symbols to the clusters is the same 
as that in Fig. 9. The inset shows the magnetization curve. In 
the configuration at m = 2/9, the interaction energy originates 
from the fourth-neighbor repulsion indicated by gray lines. 



§4. Summary 

We have been studied magnetization plateaus in the 
Shastry-Sutherland model. We have derived the fourth- 
order effective Hamiltonian, 7^g E , and then we have 
carried out partial diagonalization to get 7i^ff E for the 
subspace without first- and third-neighbor triplet dimer 
pairs, in which second-order energy takes the lowest 
value. We have checked the validity of our approach to 
SrCu2(B03)2 by using the exact diagonalization method 
for the 12-dimer cluster. It has been shown that in the 
subspace excluding first-neighbor TD pairs the corre- 
lated hopping terms, fi n _>.2n and i*2 n _>i n , lead to the 
cancellation of the second- neighbor repulsion in H*g E . 
Then the fourth-neighbor repulsion is responsible for 
the 1/4-plateau observed experimentally. At the 1/4- 
plateau, the diagonal stripe arrangement of TD's is sta- 
bilized. 
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